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ABSTRACT 


The feasibility of using radiation absorbing gaseous media is analyzed theoretically 
and experimentally verified. The radiation absorbing gaseous media like C Oo and 
NH:i tend to reduce the thermal radiation falling on the inner side of the glass cover, 
from the absorber plates by absorbing it when it passes through it. ( But these gases 
are nearly transparent for solar radiation. ) Simultaneously the gaseous medium 
emits thermal radiation uniformly in all directions, increasing the thermal radiation 
falling on the inner side of the glass cover. But the temperature of the gas being 
lower than that of the absorber plate, and the thermal absorptivity and emissivity of 
the gas being approximately equal, results in the net reduction in the total thermal 
radiation falling on the inner side of the glass cover. This reduced thermal radiation 
on the glass sheet causes lowering of the outer cover temperature, which reduces the 
convective and radiative heat losses from the outer side of the glass cover. In the 
present theoretical work, a specific type of solar collector has been designed and the 
heat collection efficiency at various outlet temperatures is calculated for collector and 
gaps with air. with gaseous media and with vacuum. For simplicity of the computations, 
walls are eissumed to be gray, highly absorbing and emitting and gaseous media to be 
isothermal, this enables to use the Hottel's emissivity charts. 
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Chapter 1 


Introduction 

1.1 Solar energy : Promises and Challenges 

• Solar radiation is incident on the earth's surface at the rate of approximately 

8x which is more than 10000 times the present world energy consumption. 

• Our energy use includes several forms of delivered energy, in particular, heating 
and cooling, light ing, mechanical motion and the industrial production. In most 
cases the useful energy is different from the input energy form. 

When considering the renewable energy sources, attention must be paid to the 
best possible match between available energy type and the desired end use energy 
type. Table 1 shows the conversion efficiency for different solar energy conversion 
methods. 

• From the Table 1 , it is evident that biomass should not be used for heating 
purpose because the conversion efficiency of the solar to biomass is very low , 
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(b) High terr.r)eraiure 
Centre: reciever 

Parabolic dish 10-25 

Pomt-focus Fresnel lens 


(iil) Biomass 


0.1-2 


and Typical Annual Average System 


Comments 

For a given system the efficiency 
decreases with temperature. 
Higher temperature collectors 
tend to cost more. 


Probably 10% is the minimum 
needed for commercial success. 
Photovoltaic systems could be 
deployed in small modules (e.g. 
on roofs). 


Shon term energy storage as 
heat. 


£nerg\^ is stored as chemical. 
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Table 2: Energy (,'onsumption Pattern of ESA in 1977 [Ul. 


Transpofiaiion 

21.2 (26.2%) 


Total 
S0.8 j 
( 100 %)^ 


i Residenitai 

I 18.2 

Residential plus i (2_.5%) 

commercial J 

29.9 ^ 

(37.0%) i 

i Commercial 
11.7 

: (14.5%) 


'8.7 (10.8%) Space heat 

2.5 (3.1%) Hot water 

1.2 (1.5%) Air conditioning 
2.4 (3.0%) Refrigerators and freezers 
1.1 (1.4%) Lights 
- 2.3 (2.8%) Other 


5.2 (6.4%) Space heat 

2.5 (3.1%) Air conditioning 
0.26 (0.3%) Hot water 

2.5 (3.1%) Lights 

1.2 (1.5%) Other 


Industnal 

29.7 

(36.8%) 




Process heat 
12.0 
(16.0%) 


- l.O (1.2%) Below 100‘C 
. 3.8 (4.7%) I00-I7rc 

2.9 (3.6%) 177-2SS-C 

2.6 (3.2%) 2S8-593-C 
1.0 (1.2%) 593-l090'C 

- 1.6 (2.0%) above 1090*C 


Other industnal 
16.8 
( 20 . 8 %) 


of the order of 1 % onl\‘. where as the conversion efficiency of solar to thermal is 


about 30% to 60%. 


1.2 Energy Consumption Pattern and Scope for 
Solar Energy Utilization 

Table 2 shows the Energy consumption pattern of USA in 1977. It is found that 
approximately 2.4 — 4.2^ energy is used for heating purpose with the temperature re- 
quirement below 100 °C. which can be easily supplied by the use of simple nontracking 
solar collectors. The energy required in the range of 100 °C to 177 °C ( about 5% of 
total energy consumption ) can be supplied by the nontracking evacuated collectors. 
Compound Parabolic Concentrator collectoi^are available that can operate efficiently 
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upto 300 °C. The total heating requirement, industrial plus residential, accounts for 
about 309c of the total energy consumption, which can be supplied by solar energy 
thermal conversion methods. 

These numbers suggest a very large potential for the application of active solar 
collectors. However, one also needs to consider the relation between the demand and 
the solar availability. The utilization of the solar equipment is best if there is year-round 
load. The role of co-generation is also to be emphasized, because when converting solar 
to electrical, most of the energy is wasted. 

Table 3 gives a guide line for the use of nonconventional energy resources in the 
coming future [29]. It lists the appropriate energy supply technology for different end 
uses. 


1.3 Survey of Solar Collector Technology 

The most important and the most expensive component of an active solar energy 
system is the collector. The survey of the sohir collectors, summarized in the tabular 
form is given below in Tables 4, 5 and 6 [14]. 

The nontracking collectors are the most economical and the best suited for the ap- 
plication of low temperature requirements. Flat plate collectors are the best developed 
collectors. A promising low-cost lightweight plastic collector has been developed by 
Andrew and Wilhelm, 1980, for application where the water quality inside the collec- 
tor is not as critical as for domestic hot water. Among the solar ponds.shallow solar 
pond and deep or salt-gradient solar pond are entirely different. The shallow water 
pond consist of a shallow horizontal water bag, insulated by one or more plastic films 
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•ble 3: A Proposed Long term Solar Energ\’ Economy [29] 





I’FHCFNIAOr Of- 

AIM’IIOI'HIAH f Nf HOY 

MAND srcion 

Ur.E ENriKlY FOMM 

ArCLICAMCJH 

TOIAL ENfJtGY USE 

SUPf^Y fCCfINOLOGY 

lidoritial and 

Low toniporalufo thorinnl 

Space healing, water heating. 

20-25 

Passive and active solar sys- 

nmercial 

energy (• 212 F) 

air conditioning 


tems, district heating systems 


intermediate temperature 



Active solar heating with con- 


thermal energy (212-572' F) 

Cooking and drying 

5 

centrating solar collectors 


Hydrogen 



Solar thermal, thermochemi- 
cal. or electrolytic generation 


Methane 



Biomass 


Electricity 

Lighting, appliances. 

- 10 

Photovoltaic, wind, solar ther- 



refrigeration 


mal. total energy systems 



Subtotal 

■ 35 


istriaj 

Intermediate temperature 

Industrial and agricultural 

7 5 

Active solar heating with flat- 


thermal energy {• .572 F) 

process heat and steam 


plnlo collectors, and tracking 
solar cojicentrators 


High tornpornturo ( 572 F) 

Industrial process heat and 


Tracking, concentrating solar 



steam 

17.5 

collector systems 


Hydrogen 



Solar thermal, thermochemi- 
cal, or electrolytic generation 


Electricity 

Co-generation, electric drive. 

10 

Solar thermal, photovoltaic, co- 



electrolytic, and 
electrochemical processes 


generation. wind systems 


Feedstocks 

Supply carbon sources to 

5 

Biomass residues and wastes 



chemical industries 


or plantations 



Subtotal 

40 


sportation 

Electricity 

Electric vetiicles. electric rail 

10-20 

Pliotovoltnic. wind and solar, 
thermal-electric 


Hydrogen 

Aircraft fuel, land and water 


Solar thermal, thormochnmi- 



transporinlinn vchir ics 


cal. or rjloctrolvlic gcMinralicjn 


Li(|uid fuels inolhnnol, 

Long distarreo land and water 

5-15 

Biomass residues and wastes 


ethanol, gasoline 

transportation vehicles 


or plantations 



Subtotal 

" 25 



100 


jRCE; Henry Kendaii and Steven Nadis, Energy Strategies: Toward a Solar Future, pp. 262-263. 
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Table 4: Xontracking Collectors 


Collector type 

.Approximate 
maximum 
operating 
temperature (®C) 

Cost (Si nr) 

Comments 

Shallow solar pond 

40-60 

160* 

(complete system, 
including storage 
for 1 day) 

Plasuc covers may need 
to be replaced every' 5 
>T or so. Needs sunny 
climate for good 
penbrmance. 

Deep solar pond (salt 

40-Q0 


Collector and longterm 

gradient I 


(includes storage) 

storage m one unit. For 
seasonal storage, depth 
should be aoout 3 m. 

Low cost, but low 
emdencv- 1 10%“20%). 

Rat plate 

(a) Conventional 
design 

40-80 

150-300"-^ 

Best knowTi and most 
developed of ail 

(b) .Made of plastic 

30-60 


coileaor T>-pes. 

(c) L'n glazed 

io-:o 

above ambient 

70-lOO"-‘' 


Nonevacuated CPC 

80-120 

150^*= 

— 

h.xed-iiit or summer- 

lo-winier tilt 




adjustment 

Evacuated tubes (with 
reflector enhancement: 
e.g.. CPCi 

100-200 

250-300^-*' 

Many opponunities for 
cost reduction oy mass 
produaion ana for 
penbrmance 
improvements through 
RSiD. 


^Cosi esnmaie in 1980. Personal commumcaiion. Solar £nerg>- Group. La^^Tencc Livermore Laborairnc 
Livermore. C.\ 9455‘j. 

^Costs of solar ponds are ver>’ siie-soecinc (e.g., sail may be tree, liner mav noi be ncedetL eic.L .A. 3-m 
deep sait-eradieni soiar pond was built m Miamisburg. OH for 35 S/m* in 1978. 

Solar P’-'^ouas Soeancaiions Guide. 1983. Published annually by Solar Aee .Mazazine. Cauicn Hill. Har- 
nsville. NH 03450. 

'^For some collectors several manufacturers are in the market, with a wide spread in quaiirv- ana pnce (not 
always correlated). Some models e.xceed the pnce range indicated here. 
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Table 5: One Axis Tracking Collectors 


Coiiecior 

Approximate 
maximum operating 
temperature (°C) 

Cost 

(S/m-) 

Comments 

Inflated cy’nncncal 
reflector 

140 

50-70“ 

Does not need continuous 
tracking, but does require 
weekly tilt adjustments: 
plastic cover may need to 
be replaced every' 5 yr or 
so. 

Parabolic trough 

300 

ISO-SOO”" 

Continuous accurate 
tracking: sensitive to din. 

Line-focus Fresnel 
reflector 

250 


May combine advantages 
of parabolic trough and of 
central receiver for 
temperatures below 

250X. 

Fixed linc-fociis 
reflector with 
tracking receiv er 

250 


Problems with din 
accumulation on 
reilectors. 


*^Cosi esiimaii ;r. 1980. Personal communicaiion. Soiar Energy Group. Lav-rence Livermore Laboratory, 
Loermore. CA ^~550. 

‘Soiar Froduc:. Spcancaiions Guide. !9S3. Published annually by Soiar Age Magazine. Church Hill. Har- 
nsviiie. NH 0345 ' 

'For some coiieciors several manufacturers are in the market, with a wide spread in Quality and pace t not 
aiwa>s correlaiec . Some models exceed me pnee range indicated here. 


Table 6: Two Axis Tracking Collectors 


Collector type 

.Approximate 
maximum operating 
temperature ("C) 

Cost (S/ 
m-) 

Comments 

Parabolic dish or point- 
focus Fresnel lens 

■ 1500 

(possibly more) 


Good if energy can be used 
directly in focal zone (e.g.. 
photovoitaics or soiar 
thermal power): otherwise, 
transponing heat to point 
of use IS problematic. 

Central receivers 

1000 

(possibly more) 

492* 

355'’ plus 
tower' 

Optical transpen of 
energy. 

FLxed-hemisphencai 
retlector. tracking 
receiver 

400 


Problems with heat 
iranspon to point of use, 
and with dirt 
accumulation on reflector. 


^Average cost of bciiosuis for Barstow solar power plant, in 1980 [Bartleson, 1981]. 
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and air layers [(’lark and Dickson, 1980]. The deep solar pond uses a thick layer of 
nonronvecting water as insulation [Tabor and Weinberger. 1980: Neilson. 1980]. Here 
coii\'ection is prevented by adding salt in such a way as to establish a concentration 
gradient, with a saltier water at the bottom. 

Evacuated and concentrating collectors are used to reduce the heat loss. While 
several attempts have been made to build evacuated flat plates, they do not seem to 
hold any promise of commercial success. The natural configuration for an evacuated 
collector is the glass tube. ,\11 of the evacuated collectors use selective coating as 
absorber because with a nonselective absorber, radiation losses would dominate at 
high temperatures and eliminating convection would not be very effective. 

In concentrating collectors, heat losses are reduced by concentrating the radiation 
incident on the aperture onto a smaller absorber. Presently, the parabolic trough is 
the favored collector for the temperature range of bdO to 300 °C. (,'ollertors with a 
line focus Fresnel lenses pose a problem for large thermal insulation. These collectors 
have been tested upto 330 .A.ttainment of temperatures above 500 °C. requires 
point-focus concentrators. The principal types are the central receivers, the parabolic 
dish and the point-focus Fresnel lens. 

(.’ollectors with high concentration reciuire accurate and continuous tracking and 
are quite sensitive to dirt accumulation, on both the reflector and the receiver. An 
interesting alternative to a concentrator with tracking reflector is a fixed reflector with 
tracking receiver. Last concentrator type is hemispherical fixed reflector with tracking 
receiver. Like other fixed reflectors it is quite sensitive to dirt accumulation and like 
the parabolic dish, it faces serious problems of heat transport from absorber to the 
point of use. 
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0 100 200 300 

Tamparitur* rO 


Figure 1: Efficiency and Temperature Range of Various Solar Collectors [*^1] 

The suitable type of solar collector can be selected on the basis of the tempera- 
ture requirement and the cost of the collector. Figure 1 shows the efficiency and the 
temperature range of different typesof existing solar collectors [31]. 



Chapter 2 


Literature Survey and Present 
Work 


In the light of the present energy crises, people have started thinking seriously in 
the direction of utilizing non-conventional energy resources and thus extensive work is 
being done to make the use of these non-conventional energy resources conventional. 
The use of the solar energy had started a long back, mainly for drying purposes. With 
the time it extended to the area of water and space heating and then photo- voltaic 
applications. 


2.1 History of Solar Energy Utilization [15J 

First instance of the use of solar energy is found in a newspaper of France ( Seventeenth 
century ). It states that .A.rchimedes set fire to Marcellus’ navy by means of a burning 
glass composed of small square mirrors moving every way upon hinges which when 
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placed in the sun’s rays, directed them upon the Roman fleet so as to reduce it to 
ashes at the distance of a bowshot’ [15]. 

For the next 1000 years, most of the devices invented to use solar energy were 
constructed for amusement purpose. Solar Fountain of de Caus, used several glass 
lenses to concentrate the sun’s rays on an airtight metal chamber partially filled with 
water. The sun's rays heated the air, thus forcing water out in a small fountain. 

The French Scientists George Buffon constructed the first multiple mirror solar 
furnace in 1747. Lavoisier, was able to melt platinum at 1780 °C, with a special 130 
cm. diameter lens filled with alcohal, in the solar furnace. 

Pioneering work with flat plate collector was done during the middle of the eigh- 
teenth century by the Swiss scientist Nicholes de Saussure. He designed a solar oven 
consisting of glass plates spaced above a blackened surface enclosed by an insulated 
box. 

A concentrating type of solar cooker was described in an article in Scientific Amer- 
ican in 1878 by W. Adams of Bombay, India. Solar cooker that employed parabolic 
mirror to reach high temperatures and used a means of heat storage enabling food to 
be cooked after sundown, was developed by Dr. Charles Abbott. 

Development of solar powered engines was begun in the latter part of nineteenth 
century by Augustin Mouchot in France and .John Erricsson in the United States. A 
water pump, energized by a large solar reflector, was built by A. G. Eneas in 1901. In 
1908, H. G. Willsie and John Boyle, Jr. built a 15kW solar engine, but the project was 
financial failure. Frank Shuman made solar powered engine using flat plate collector 
which was commercially successful. A solar still was built by J. Harding and C. Wilson 
for fresh water production in desert areas. In 1950, Chaplin. Fuilen and Pearson 
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developed the first solar cell. In 1959. the first successful Vanguard satellite carried 
solar chips providing power for space applications. 


2.2 Literature Survey 

Though, we had different type of solar collectors by the end of the ninenteenth century, 
they were not used in day today life because of many problems associated with them. 
They were costly in comparison to conventional methods as well as their use was quite 
troublesome. 

The main problem with the solar energy utilization is its low intensity and uncer- 
tainty regarding availability. Many people have worked in the field of energy storage 
devices to overcome the difficulty of the uncertainty. Due to low intensity of the so- 
lar energy, heat collection efficiency decreases rapidly with the outlet temperature. 

Use of multiple covers to reduce the front heat losses has been analyzed by many re- 

a 

searchers. Use of concentrating reflectors is costly affair, but is important in the case 
of higher temperature requirements. According to Collares Pereira and Sequeira, 1982 . 
Compound Parabolic Concentrator reflector with low concentration and non-evacuated 
receivers may be economical for temperatures around 100 °C [14]. O'Gallagher et al 
1982, achieved efficiency of about 50% around 200 °C with fixed tilt evacuated col- 
lector using glass tubes shaped as compound parabolic concentrator [14]. The next 
possible step was to get the optimum orientation of the collector, so that maximum 
energy can be harnessed. Many researchers have done work in this direction. [33] The 
next step was to reduce the heat losses from the cover. People tried to use the selec- 
tive coating on the absorber surface so that it absorbs maximum of the solar radiation 
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falling on it but emits the least thermal radiation back to the cover. It was found 
difficult to use these selective coatings in everyday life because selective coatings need 
very careful operation to protect them from damage. The other way to reduce loss is ro 
reduce the convective loss inside the cover gaps, which will result in increased collection 
efficiency because of the increased resistance to the heat flow tow^ards the atmosphere 
through the cover. For this purpose people are trying to evacuate the gap between the 
two cover sheets [28], [32]. Glass is the most popular cover material used because of irs 
transmission properties, but it can not withstand any compression as it is very brittle. 
Thus evacuating the cover gap is not an easy task. The natural configuration for this 
is the tube form. .A.11 of the evacuated collectors must use a selective absorber coating 
to reduce the radiative losses at high temperatures. People are using different arrange- 
ments to make it possible but no one has got very satisfactory results. The third way 
to reduce the losses may be the use of the radiation absorbing media either fluid or the 
liquid. The present work is to check feasibility of the use of the gaseous media for the 
purpose of increasing the heat collection efficiency along with the temperature. The 
use of this property of the radiation absorbing gaseous media has been studied in a 
different way and also in a different application Tn solar thruster’ by Venkateswaran 
et al [5]. The emissivity of the gases like CO^-, N Hz etc. has been theoretically and 
experimentally obtained by researchers like Hottel [8]. Methods have been evolved for 
calculating radiant heat exchange in a gas filled enclosure with gray walls and isother- 
mal gas by Hottel and Mangelsdorf [6] and Eckert. Bevan and Dunkle [3] had given 
the method for calculating radiation exchange in an generalized ceise using the band 
measurement data.s. Mea.surements by Howard, Burch and Wiliams [1] give the band 
absorption for the complete spectrum of the gas at room temperature. Hottel and 
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Cohen [7] have worked in the direction to get the allowance for nonuniformity of the 
gas temperature. Edward [1] also suggested the methodology to solve the case when 
enclosure is of nongray walls and gas is isothermal. 


2.3 Present work 


The present work is aimed at the study of the effect of the use of the radiation absorbing 
gaseous media in flat plate solar collectors. A computer program has been made for 
the purpose of the theoretical study. It is made for tracking collector in polar mode 
and for nontracking collector. The geometry of the collector is so designed to extract 
the maximum amount of the energy. ( For details refer Chapter : Experimental set-up 
) Emissivity of the gaseous media have been taken from Hottel’s Emissivity Charts. 
Program has been made for transient study of the collector. A model of the collector 
had been fabricated and tested for Ammonia. The analysis of the effect of evacuating 
the whole collector has also been studied and the gain in the temperature and the 
collection efficiency has been compared for different cases. 

The results obtained show sufficient gain in temperature and collection efficiency at 
higher temperatures in the case of absorbing media. At large flow rates, there is negligible 
gain in the efficiency and outlet temperatures. The gain is increasing with the increase 
in the working fluid outlet temperature. The experimental I'esults also showed increase 
in the efficiency when gas (NHs) was used. Thus the computational results are verified 
qualitatively. Since the use of these absorbing media is neither a costly affair nor it causes 
any severe problems except leakage of the gas, it can be implemented easily in all types of 
the solar collectors where the leakage can be avoided. 



Chapter 3 


Fundamental Concepts 


3.1 Solar Energy 


3.1.1 The Sun and Solar Energy 


The source of solar energy is the Sun, which is a gaseous star (approximately 80% 
Hydrogen and 18% Helium ).For all practical purposes, it may be considered to be ra- 
diating energy at an effective temperature of 5760 °K. The spectrum of the wavelength 
of radiation stretches from 0.29 to 4.75^m. The maximum radiation intensity is found 
at wavelength of 0.50;um. In order to calculate the beam and diffuse solar radiation in- 
tensities on different days and at different times of the day at different locations on the 
earth, the solar geometry and the location of the place of interest need to be defined. 
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Figure 2: Solar Altitude and Azimuth angles with reference to local horizontal plane 
for North Hemisphere. [l2] 

3.1.2 Solar Geometry 

It is adopted that the sun moves on an imaginary spherical surface called celestial 
sphere and this sphere is sufficiently large so that the sun can be considered as a point. 
The celestial sphere is assumed to he centered at the site of interest on the earth's 
surface. .Since the sun is moving on a sphere, (i.e. two degrees of freedom ). two angles 
are recjuired to define the instantaneous position of the sun. 

[1: Solar .Altitude angle a 

It is measured from the local horizontal plane upwards to the center of the sun. It 
is zero at sunrise and sunset and maximum at solar noon. (See Figure 2 ) 

[2l Solar Azimuth angle 

It is measured from in the horizontal plane between the north-south line and the 
projection of the sun's rays onto the horizontal plane. It is zero at the local solar noon, 
positive before noon and negative after noon. (See Figure 2 ) 
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The solar altitude and azimuth angles are related to the solar hour angle the latitude 
L. and the solar declination Ss- The solar altitude angle a. is given by 


sin{a) = sin{L) .sin{6s) + cos{L) cos{5s) cos{Hs) {sin{(y) >= 0) (1) 

The solar azimuth angle is given by 


sin{as) = 


co.s{6s) sin{Hs) 
cos{a) 


Other angles which come into calculations are. 


( 2 ) 


[1] Latitude L 

It is the angle between the earth's equatorial plane and a line from the center of 
the earth to the site of interest. It is positive in the north of the equator and negative 
in the south. 


[2] Hour angle Hs 

It is given by 15 ° times the number of the hours from local solar noon. It is zero at 
noon, positive before noon and negative after noon. [ The signs of the angles are based 
on the usual positive angular displacement rule in right hand coordinate systems.] The 
expression for Hour angle is. 


Hs = 12.0 — corrected standard time — ETjQO (3) 

'Where ET is the Equation of time, which is the difference between standard clock time 
and the solar time. ( This happens because the length of the solar day varies due to the 
tilt of the earth’s axis, the nonspherical shape of the earth and the orbital eccentricity. 
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) It is expressed in minutes and its values are tabulated in Appendix A, Table 7 
An empirical equation for calculating the ET values is given by DufBe and Beckman 
1980 [14], 


= 9.87 .sm(2B) - 7.53 co.s(B) - l.o sin{B) (4) 

Where B is and Nday is the number of the day i.e. 1 for .January. 

1 and 365 for december, 31. 

The corrected standard time for west of Greenwich, is defined as. 


Corrected standard time = Standard clock time 

+ {time zone meridian — local site longitude )/15 


[3] Solar Declination Ss 

It is the angle formed by the line from center of the earth to the center of the sun 
at noon on a particular day and the plane containing the earth’s equator.(See Figure 
3 ) Cooper suggests the approximate equation 


36o 


V 


J 


( 5 ) 


For engineering purposes, the declination is cissumed to be constant over a day. The 
declination is considered positive when sun is in the northern latitudes and negative 
when in the southern latitudes. On .June 21, summer solistice, sun is at its northern 
extreme position hence declination is maximum positive and is the longest day in 
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Collector site at latimde o 




Figure 3; Solar Declination angle 


northern hemisphere. On December 21. winter solistice. sun is at its southern extreme, 
hence declination is maximum negative and is the shortest day in northern hemisphere. 


3.1.3 Solar Incidence angle di 

The solar incidence angle is the angle between the normal to the irradiated surface and 
a line collinear with the sun’s rays. The solar incidence angle for any fixed surface is 
given by 


co.s{di) = sin{6s)[sin{L) co.s{/3) — co.s(L) szn(/3) co.s(ac)] (6) 

+co.s((55) cos(Hs) [cos{L) co.s{0) + sin{L) sin{fi) cos{ac)] 

+co${6s) sin{3) sin^a^) sin(Hs) 

where 0 is the surface tilt angle from local horizontal and is the surface azimuth 
angle defined analogously to the. solar azimuth angle a^. 

The angle 0 is positive for a south facing surfaces and negative for north facing 
surfaces. The negative value of cos{0i) indicates that 9i > 90 ° and the sun’s rays do 
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not strike the surface. 

3.1.4 Incidence angles for Tracking Collectors 

Tracking is required, mainly in the case of concentrating type of collectors, to maintain 
the sun within their field of view and to minimize the incidence angle. Either one or 
two degree of freedom modes can be used according to the requirements. 

For perfect tracking, co.s(0, )=l.O Single axis tracking can involve either a horizontal, 
east-west axis or a horizontal, north-south axis of rotation. Also a polar mount can be 
used which employs a north-south axis tilted up at an angle equal to the local latitude. 
The incidence angle for these modes is, 

1 East- West axis Horizontal 

co.s(^i) = (1 — cos^{a) sin^ (7) 

2 North-South axis Horizontal 
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co.s{di) = (1 — co . 5 !"( q ) sin^{as))^^^ ( 8 ) 

3 Polar axis 

CO.S(0;) = COsiSs) (9) 

3.1.5 The Extraterrestrial Solar Radiation 

The mean annual value, of the radiation intensity normal to sun’s rays incident upon a 
plane surface situated in the outer limits of the earth's atmosphere is called the solar 
constant. It varies inversely with the square of the distance from the center of the 
sun. The earth axis of rotation is tilted 23.5 ° with respect to its orbit around the 
sun. therefore, solar constant varies with location ( latitude ). The present accepted 
standard values of the solar constant [determined by Thekackara and Drummond ] is 
1353 W/m". The spectral distribution of the extraterrestrial radiation. Figure 5 is 
Ultra-violet region range 0.00 to 0.38;xm 7.0% 

Visible region range 0.38 to 0.78/im 47.39c 

Infra-red region range 0.78 to oofim 45.7% 

The value of the extraterrestrial solar radiation changes with the position of the 
sun because of the elliptic orbit of the sun. ( Appendix A.Table 8 ) The equation 

for the value of the extraterrestrial solar radiation is 

h = he [l.O + 0.034 cos (10) 

The value of the extraterrestrial horizontal solar intensity at any time is given by 
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Figure 5; The Extraterrestrial Solar Spectrum and Molecular Absorption Bands of Gases 
in Atmosphere 


loh = h [•sm((5s) sin{L) + cos{5s) cos{L) cos{Hs)] (11) 

A. 

3.1.6 The Terrestrial Solar Radition 

As the solar energy filters through the earth’s atmosphere, it heats the atmosphere 
by molecular absorption and a part of it is scattered in new directions by molecules, 
water droplets, ice crystals, dust particles etc. Loss of energy through scattering and 
molecular absorption occurs mainly in visible range, whereas, loss in the infrared range 
is caused by absorption by water vapor and CO 2 present in the atmosphere (Figure 
5). 

This results in a net loss of 0 to 80% of the extra-terrestrial solar heat flux. The 
remaining, 0 to 80%, arriving at earth’s surface, now, consists of two components, first 
is direct or beam radiation and the second is diffuse radiation. The direct or beam 
component appears to come directly from the sun and diffuse component is basically 
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from scattered and re-emitted energy by molecules in the atmosphere. The amount 
of the energy incident upon a terrestrial surface, pointed directly to the sun and its 
distribution in direct and diffuse components, depends upon atmospheric conditions. 
The terrestrial beam solar radiation intensity is given by. 

h = Io[(i{F't) Aj -r &(Ar)] (12) 

Where Kj is the Clearness index or atmospheric transmittance and a(Kj) and 
h{Kj) are empirical coefficients (.Appendix .A.Table 9 ). 

The Terrestrial horizontal solar intensity at any time is given by. 

h = I oh A'r (13) 

The Terrestrial horizontal diffuse solar intensity at any time is given by, 

= Ik - h cosill) (14) 

The diffuse solar radiation intensity on collector surface at any time is given by, 

he = hh co.s^(|) (15) 

3.1.7 Sky Radiation 

Infrared radiation exchanges occurs between solar collector and the sky. Since the 
temperature of the deep space is about 2.7K, the effective sky temperature for radiation 
is usually below that of the ambient air temperature Ta- The net exchange Q surface-sky 
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between a surface ( gray ) on the earth at temperature Tg . and emittance and the 
sky at eifective temperature Tgky can be written as 

Qsurface-sky = eg<r(T^ - Tg\^} (16) 

The sky radiation depends on the radiative properties of the sky (e.g. humidity 
and cloud cover ) and many other secondary phenomena. An equation for Qsurface-sky 
including principal climatic effects is given by deJong. 

Qsurface-sky = (0.39 - 0.009e) [l - UsiCC )] (17) 

+4tsCrT^{Ts — Ta) 

Where e is the vapor pressure of atmospheric water vapor in mm Hg and CC is the 
cloud cover expressed in tenths . It is 1.0 for full overcast day and 0.0 for cloudless 
day. a 3 is the empirical constant (Appendix A.Table 10 ). 

Thus the infrared radiation coming from the sky is, 

Qti Q surf ace— sky 

3.2 Transmissivity, absorptivity and reflectivity of 
glass sheet 

The transmissivity absorptivity Og and reflectivity Pg of a glass sheet depends upon 
the refractive index n, absorption coefficient 7 ( the optical properties of glass) and 
the sheet thickness Ag. 



(18) 
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Figure 6: Spectral Absorption Coefficient for Different Types of Glasses [b] 

3.2.1 Absorption coefficient 7 

For microscopically homogeneous transparent media. Beer's law states that 

Ix{x) = (19) 

where x is the layer thickness and 7,\ is the absorption coefficient at wave-length A. 
Thus, for a single pass through the sheet, transmittance r and reflectance p is 


p'j. 

p'w 


— 


g COS{ ) 


.sin^(di — 

Or) 

sin^(0i -4- 

Or) 

tarp{di — 

Or) 

taTp{di 4- 

■Or) 

{n-lf 



(n + 1)'^ 


(for 9i > 0.0) 
(for Oi > 0.0) 
(for 6i = 0.0) 


( 20 ) 
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Figure 7: Ray Trace Diagram for Glass Sheet 

In the case of sheet there is multiple reflections from the two interfaces, hence each 
of the. Tg. ag and pg is expressed in the form of an infinite geometric series as a function 
of T . a and p . (See Figure. 7 ) These equations are known as Stoke’s equation. In the 
case of high absorption coefficient, the refractive index in the expression of reflectivity 
becomes complex number, i.e. 


rieff = n + 


Xkx 

drr 


( 21 ) 


It is very difficult to account this elfect in exact manner, as the data for kx is not 
available. (See Figure 6 ) But some approximations are required for the reflection 
of the thermal radiation where absorption coefficient becomes very high. This effect 
can be incorporated by increasing the value of refractive index by some appropriate 
amount for calculation of reflection from the surface. 


Tg 


.'2 


'^2 


1 + 


(1-p) 
1 - 


( 22 ) 
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When both component of polarization are equal, the average reflectance and average 
transmittance are 


0.5[p(p'^.r^) + plpjp q|)] 

(23) 

0.5[r(p'_L,rl) + r(/9||,r||)] 

(24) 


Once the values of effective transmittance and reflectance are calculated, the total 
absorptance can be calculated as 


- pg ( 25 ) 

3.3 The Emission of Radiation by Glass 

In transparent materials, the emission and absorption of radiation is a bulk phenomena. 
Interaction of the simultaneous emission and absorption of radiation throughout the 
volume determines the rate at which radiation appears to be emitted by the surface. For 
a nonscattering dielectric material this rate depends upon the temperature, thickness 
of sheet, refractive index and absorption coefficient of the material. 
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3.3.1 Volume Emissive Power 

Arrording to McMahon [8], the spectral volume emissive power is a measure of the 
monochromatic power radiated at a given temperature by a unit volume of material 
into unit solid angle of space. 

For isotropic materials, radiation is a diffuse process which proceeds uniformly in 
all directions. 


. IxnnVsx 



/t 

For a gray material, total volume emissive power is 

. inHVB 

j ^ — _ — 

TT 


(26) 


(27) 


3.4 Geometric Mean Beam Length 

In the case of absorbing and emitting media within an enclosure, the radiation from 
any differential surface element travels different path lengths to reach the other surface. 
The effect of these different path length can be expressed by a single mean path length 
known as Geometric Mean Beam Length. It is only dependent on geometry of the 
surfaces. Mathematically, It is defined as, 

rm _ 1 f f cos{l3j) cos{^k) dAk dAj 

- AjTj Ja, Ja. 

Where is the path length and /? is the angle made by D’ from the normal to the 
surface element dA. 



CHAPTER 3. FUNDAMENTAL CONCEPTS 


29 


The value of the geometric, mean beam lengths for two parallel plates of side a and 
b separated by distance c and for two perpendicular plates of sides (a and b) and (6 and 
c) have been calculated by Dunkle and Hottel [20]. The values calculated by Edward 
are given in Appendix B, Table 11 , 12 


3.5 Gas Emissivity and Absorptivity 

Until the early 1920s gas radiation was treated as a correction to the convective flux 
from a gas. The problem of evaluating the radiant heat transmission from gas flames 
received attention because of the practical importance in the design of the industrial 
furnaces and internal combustion engines. 

A change in the occupancy of a particular energy level of the atom or molecules of a 
gas. can be viewed as emission or absorption of radiation in a narrow wavelength region. 
This process is known as line emission or absorption. This type of selective radiation 
spectrum is the characterstic of the gases. For monoatomic and diatomic gases, the 
capacity to emit and absorb radiation is insignificant, but in the case of poly-atomic, 
gases like CO 2 . water vapor. SO 2 and NH:i, they have considerable emissive power and 
high absorption capacity. 

Gases with polar molecular form i.e. asymmetric molecules, exhibit energy trans- 
formations in the wavelengths of thermal radiation. The typical thermal energy trans- 
formation occurs as a result of changes in the vibrational frequency and rotation. This 
is manifested in a strong peak of energy emitted at a wavelength corresponding to the 
vibrational transformation, with multiple rotational energy emissions at wavelengths 
slightly different than the wavelength for vibrational energy. For this reason energy 
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emitted is usual!}’ described as a band of energy over a group of wavelengths from a 
rotational-vibrational transformation. 

The result of this type of energy interaction is that gaseous absorption or emission 
occurs in discrete bands, giving a highly nongray behavior to the gases. 

3.5.1 Gas Emissivity 

The term gas emissivity has a meaning only in reference to emission from an isothermal 
gas shape to a specified portion of its bounding surface. It is the ratio of the incidence, 
on that surface portion, from gas lying within a specified field of view to the incidence 
from a black emitter at the gas temperature over the same field of view; and it is 
plainly dependent on the shape of the gas mass. The standard emissivity of a gas is 
that corresponds to a fixed path length. 

Similarly, the standard gas absorptivity is the fractional absorption, by an isother- 
mal gas. of blackbody radiation emitted from a bounding surface element through a 
fixed path length. 

Similar to Kirchhoff's law for surfaces, 

^ ^ ^ _ g-(.-ypL"*) ( 29 ) 

The total emissivity and absorptivity are, 
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^mix 


r W- a.,^A 


— f 

aT^ Jo 


a,s<^A 


Where Tg and Tg are the surface and the gas temperatures respectively. 


(31) 


3.5.2 Interrelation between Absorptivity and Emissivity 

In terms of band black widths, the definitions of gas emissivity and gas absorptivity- 
are 


,n.^^Tg) = 


(32) 

a”^^^{Tg,Tg) = 


(33) 



(34) 


Where Ag, is the equivalant black width of a band. 

Similarity between the above two equations, suggests the possibility of expressing 
absorptivities in terms of emissivities. If the emissivity is evaluated at the temperature 
of the surface, the right hand side of the Equation 33 and Equation 34 becomes 
identical with the exception that Ab, is evaluated at Tg in the former and at Tg in the 
latter. Assumption of a constant absorption coefficient within a band gives 

Ab. = (35) 

Using the temperature dependence of Au and a 
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a{T,.T,,pL^) 



(36) 


3.5.3 Gas Emissivity Charts 

In view of the complexity of the theory and the associated uncertainty in the calculated 
values of gas emissivities. total radiation measurements are assumed to be most reliable 
and convenient source of information. 

In spite of the highly nongray behavior of the gases, gray gas approximations are 
applicable in the case of energy exchange in systems where the walls are either black 
or nearly black. 

The variables required to define total gas emissivity and absorptivity are 


e-- = /i(L-,p,P.C,r,) (37) 

= f2{L^,p.P,C,TgMs) 

(38) 


Where C is the composition of the remainder gas. In most cases effect of C is negligible. 
If there is more than one emitting species in a gas. a correction factor is required to 
account for the overlap among the emission bands of the different molecules. 


3.5.4 Emissivity Chart for CO 2 

Radiation from CO 2 principally comes from bands at about 2.64 to 2.84. 4.13 to 4.50 
and 13.0 to 17.0 microns. Rubens and Landenburg 1905 , Hottel and Mangelsdorf 
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Figure 8: Emissivity Chart for CO 2 at 1 atmospheric total pressure [ 8|. 


1935 , Hottel and Smith 1935 and Eckert 1937 made total radiation measurements 
for C O 2 either in pure form or in mixture form. Emissivity chart for CO 2 prepared by 
Hottel is shown in Figure 8 . 

The absorptivity and emissivity are correlated by 


(39) 

3.5.5 Emissivity Chart for NH^ 

Total radiation from ammonia vapor-nitrogen mixture has been measured by Port. 
The results are presented in the form of chart. The emissivity chart for NH^ prepared 
by Port is shown in Figure 9 . 

The absorptivity and emissivity are correlated by 
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Figure 9: Emissivity Chart for A' H 3 at 1 atmospheric total pressure [ 8J. 


/ 'T' \ / 'T \ 

f (401 


3.5.6 Equation of Energy Transfer in Nongray Isothermal 
Gas 

For analysis of the energy exchange in an enclosure that contains a nongray gets, the 
frequency band of interest is divided into absorption regions and window regions. In 
window regions, the gas is totally transparent. This causes two different types of 
expressions for energy transfer: one for regions of absorptance and other for the window 
regions (See Figure 10 ). 

If the enclosure consists of n surfaces, each of which is assumed to have uniform 
radiousity, the equation of energy transfer for each wall for absorption region is. 




CHAPTER 3. FUNDAMENTAL CONCEPTS 


3o 



WkveiengTn 


Figure 10: Band Division for Energy Exchange Calculations 


Gu.jAj = ^ ^ ^bu,{Tg)FkjAkO.;\i^^ (41) 

k=i k=i 

for j =1,2,3 n 

for window region is. 


G^jAj = ^ JuikPkjAkr^^ (42) 

<:=1 

for j =1,2.3 n 

Where Jj is the radiosity of the jth surface, defined as 

— PjF^j + (-^ 3 ) 

These sets of equations must be solved for all absorption and window regions and 
then results are summed up. Thus it is a very exhaustive technique, but unfortunately 
this is the only technique when walls are reflective. When the gray gas assumption are 
valid, the above equations are used with all quantities as total quantities. 


GiA^ = f: + y W,{T,)F,iAtc 


(44) 



Chapter 4 


Mathematical Modelling 


4.1 Assumptions 

[1] Walls are gray and their emissivity and absorptivity are near to unity. 

[2] Gas volume is isothermal. 

These assumptions justify the use of Hottel’s Emissivity charts. 

[3] Pressure of the gas inside the gaps and box remains constant. 

[4] No scattering takes place either by gas or glass volume. This makes 
the extinction coefficient equal to the absorption coefficient. 

[5] Glass properties (absorption coefficient and refractive index ) are 
wavelength independent, since the wavelengths are grouped in three 
regions, namely thermal, infra-red and solar. 

[6] Beam solar radiation becomes diffuse when reflected by glass surface 
or walls. 

[7] Since absorber plate is made of copper sheet, it’s temperature is 
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uniform throughout the plate. 

[8] The heat transfer to the working fluid takes place only from the one 
side of the channel, which is formed by the absorber plate itself. 

[9] Reflectivity of the ground is negligible. 

4.2 Notation for walls 

The walls are numbered as. 

1 ; For upper cover 

2.3 : For walls never facing sun’s rays. i.e. walls in the direction 

perpendicular to surface azimuth line in the case of tracking collector 
and facing south in the case of nontracking collectors. 

4 : For wall getting solar radiation after solar noon, also denoted by 1. 

5 : For wall getting solar radiation before solar noon, also denoted by r. 

6 : For lower plate, also denoted by b. 

Covers are numbered as 1 for top-most and Ncover for lower most. 

Strips are numbered as 1 for lower most and iVi for upper most. 

4.3 Solar Irradiation on Different Walls 

Solar radiation falling on the collector cover consists of three components, namely beam 
or direct, diffuse and radiation reflected from ground. Intensity of beam solar radiation 
If, is given by Equation 12 . The intensity of the diffuse solar radiation is given by 
Equation 15 . The total amount of diffuse solar radiation on collector cover is. 
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Qin = Idc A^ (45) 

4.3.1 Direct Solar Radiation on Different Walls 

Distribution of direct solar radiation on different walls depends on the angle of incidence 
( 9i^ measured from the horizontal plane ) of the rays. The distribution have been 
calculated for 6^ = 0 to | as it remains same for = I to except that the left and 
right wall irradiations are interchanged. 

If \dh < f then direct solar irradiation on left walls is always zero. 

For ig = 1,2, ... , N^^ver - 1 , = 0.0 

If [Of, < tan-^ (^^^7^)] then 

Qt,, = hr:\s7.nih)A, (46) 

Qt,. = hT:^,stn{h)A6 
Qfn„. - 0.0 
= 0.0 

for j = 2,3, ,( Meeker - 1) 

eS ,, = 0.0 


_1 f igAq 


> Oh > tan~ 


tg-l)A„ 

La 


, for ig 


N 

.... , X V J-Q 
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For j = 1. 2. 


- 1): 


Qt,., = hr:^, {sin{9iMb-ij-l){\ap + S,)L,cosie,}^ 

= h ^,.,L,cos{h) 


For j = ig. 


= hrl^j {sin{dh)Ah - [ig - \ ){lAgapA Ag)Li,cos{9h)) 

= h <j>i {sin{eh)Ab - {ig — l){A^gap + Ng)LbCos{6k)) 


For j = (ig+l), (ig + 2). 


V + 1 

t ’ cover 1 ^ -5 


Qt, = 0.0 


Qfns.. = 0.0 


0.0 
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Where is the effective transmissivity ( without taking reflections into account but 
absorption by the gas in the covers is considered ) of (j — 1) number of glass sheets 
placed in parallel. It is defined as. 


'e.l 


= 1.0 


(ol) 


^sb 

/ « ■> — 


^sb 

‘g 


'e,3 


so ^SD ^77ltX' 

'g ‘e.2 '1,16 


for j = 3. 4. 


X 


+ 1 

rft = 


^sb ^sb 

‘g ^’-2,16 


If 


\tan' 




ver ^ naybrj ^ cov ^ r )^fj-rLi: 


La 


^)>ef,> fan-^ 


^ cover t)(.l\<3Qp-j~^g) 


La 


■)1 


then. 


For ig — 1.2. 


QiC = hr:% {sin{9,Mk-{ig-l){A,,p + A, )L, cosier)) (o2) 
~ L X.ig+l Agn-pLl,COs{9h.) 


Qin., = 0-0 (53) 

Qinsr ~ '^e%cover^l ( .SZtl ( } Af, — {Ncovev “ ^){^gap + A^)L^€OS(0/i)) 
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If [0/i > tan~'^ I' , ^^cov.r j j 
For =1.2 N^over- 

= Ib isin($/AAf, - (ig - l)(Agap A A^)L(,cos(9k)) (o4) 

^intg,r ~ ~e,ig+l ■-4gap-I'bCOs{d h) 


Qins.b ~ ”e.Xr-..r5r+l — ((A'cover — I)^gap A Acover-^g A L^) L^CO.s(0^^p') 

QtllB.r ~ ^e..\%,M,:r+l Ll,LcCOS{9h) 


4.3.2 Total Incident Solar Flux on Different Walls 

The beam or direct solar radiation falls only on three sides of the box and gaps, but 
absorptivity being less than unity, causes the reflection of the radiation resulting in the 
solar irradiation on other walls also. Since even the direct solar radiation changes its 
characterstic to approximately diffuse, the reflected solar radiation have been assumed 
to be uniformly diffuse in all directions. 

For the irradiation of the solar radiation on the walls of the box, equations are 


G SO 

BA 


G SO 

Ba 


E(1 - rstrrct^ 

j=2 

6 

\ IP ^so I ^so\ p mix^ 

2^0 - ‘Bag ^B.J + Pg )^12 7-5,21 

j=3 


G SO I /~i' so rp ^rr.ix^° 

B.t A 12 'S.12 


( 56 ) 
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/I so \ rp j^^so . ^so \ jp r^so , so rp 

‘B3j ^N.j'^Pg '5,31 ^S,l ^ ^S,T^13 '5,13 

j=2,j^3 

6 

/'^so \ ^ /I -SO \ rp ^mix'^ r^so , ^so\ rp ^Tnix'^ /^^so , jp ^mix^'^ , p- 

''-'^BA — 2-/ A — CtB,j)^ j-i ' BA: ^-^B.J A Pg ‘BAl ^ BA N U b.T ^ ‘ BA4 k 

J=2.Jr^4 

6 

\' ^ /I ^so \ rp ^rnix^^ /^so , ^so\ rp ^mix^^ f^^so * so rp ^ 

B.5 — 2 -/ A — OcBJI^jS ‘ B.5j ^B,j ■'■ Pj )-ri5 'B,51 ^B,l ■*■ ' B, 13 O 

i=2.j¥5 


CA 


h 


G so V^/ 1 ^ -SO \ rp ^niix'*^ r^so . ^so\ rp /^so , so rp ^mtx'^ s r^sh 

BS - 2L.U “ '5,6J ^5,J + Pj )^16 '5.61 ^5.1 “ ^-^5.7^16 ' 5,16 ^ Q :r. 


^^B.b 


Where GJj- is the diffuse solar radiation coining to the box through the upper side of 
the glass sheet. 

For single glass cover. 


Gb°t = Q 


sci ^sd 
<fin ‘ g 


(57) 


For multiple glass cover. 


/^^'so sd 

^ 5,7 - ^^Ncov.r-lS^g 


(58) 


For the irradiation of the solar energj' on the walls of the gap. equations are 


6 


j^SO 

XX 

E(1 - ct/' GS., + AN (G& + pM.J 

j=2 

(59) 

/^SO 

= 

6 

"V^/l ^ “SO \ rp ^rnix^^ j^so i ^so\ jp ^mix^° /^so 

/ > ( ^ ^ig,j) ^gJ'^ ^ig,2j ^igj Pg ^ig,21 ^ig,l 

j=3 



+ 

f.,..2 an (g'4!l + Gfe + p,Qt.J 


/^SO 

= 

6 

V~^ /I ^ -50 \ p _7nzx^^ j^so 1 ^so\ rp mix^^ r^so 

Z_j ^ig,3j ^igjj ^ Pg ^ig,l 

j=2,j?^3 
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+ 

rp ^TTiix' 

° (Gii + G 
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P3QtL 
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II 

6 

Ed- 


igAj 
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p7 

\ IP ^mix‘^ 

)rzgP4 ‘^gAl 

^igA 


j=2,j7^4 
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^mix^^ f rp 
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^rmx ' ' 
‘ igSj 
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.so 

‘3.T 


Where G\g°j and are the diffuse solar radiation coming to the gap through 

the upper and the lower glass sheets as a result of transmission of the solar irradiation 
on the lower and the upper sides of the respective glass sheets. 


4.4 Thermal Irradiation on Different Walls 

Every wall emits thermal radiation at all temperatures and simultaneously absorbs 
the thermal radiation. The extent of emission and absorption is determined by the 
properties known as emissivity and absorptivity. In an enclosure walls radiate, absorb 
and reflect the radiation. Thus the irradiation on any wall consists of multi reflected 
and direct radiation from other walls. 

For the box. irradiation of thermal radiation on different walls can be calculated as 


^th 


Ed- 


J=2 


a 


AFb, 




‘B,jl 




( 60 ) 
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-L 

+ 

+ 

+ 

+ 

Where j is the diffuse thermal radiation coming to the box through the upper side 
of the glass sheet. 

For single glass cover. 


P ( tkpHh r-^'th , p)l \ 
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T ^mtx 
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E mix rr 
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th 
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(61 -I 


For multiple glass cover, 


^^B.T — ^Ncot.,r-1.6'5 


(62) 


For the gs-p, irradiation of thermal radiation on different walls and glass surfaces 
can be calculated as 


pith 


j=2 


+ (tT- 

i=2 




(63) 



CHAPTER 4. MATHEMATICAL MODELLIXG 


46 


+ J 
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Gg,5 = 

J=2 

J='2 

+ (pf Gt, + G';> 4- g-:.^ ) 

i" ^ig^SS \Pg ^^ig.6 ^^ig.L 

6 

t rpmix \ mix jt 

^ig 2— ^ 

j=l J^5 

g“, 6 = i;(i-c;v)Fi,,s-”5'‘6t^ 


J=-2 




j=2 


I ^ig,l6 ‘ ig,l6 \Pg ^igA ‘ ^ r^g 

5 

( rpmix \ ^ inix rp 

ig / ^ ^ia.6i 
J = 1 


Where and GG^^^ are the diffuse thermal radiation coming to the gap through 

the upper and the lower glass sheets as a result of transmission of the solar irradiation 
on the lower and the upper sides of the respective glass sheets. 


4.5 Heat Transfer to Fluid in Flow channels 

Rectangular channels are attached to the absorber plates, hence heat is transferred 
from the absorber plate to heat extracting fluid by the side, which is formed by the 
absorber plate. The channels are made up of G.I. sheet, hence the heat transfer by 
remaining three walls of the channels can be neglected {as conductivity of G.I. sheet is 
not good to increase the temperature of these walls). 
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When the fluid flows inside the channel, it is heated up as it moves forward. Here 
analysis for one channel is given, which can be used for all channels by putting appropri- 
ate channel length 1 and fluid inlet temperature 14,^ . The absorber plate temperature 
of the face receiving solar radiation Tp^ is assumed to be uniform through out the plate 
and temperature of the inner side Tp^ is assumed to be uniform for each channel indi- 
vidually. 

The energy balance for a differential element of length Ox at a distance x from the inlet 
leads to the ec^uation 

h Lr_hdx(Tp, 

or, 

^r/(r) 

(Tp, — ) 

On Integration, this gives 

Tp, — Tj(x] 

Applying the inlet condition results 

Tu 

The value of the Tp, can be calculated as. 

Rearrangement of the terms yields, 

T 


Pj Tj^ P/, 


rp ^pPj^^chyjCpf ^ ^ 

iiyF-- 



( 65 ) 
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These equations are solved iteratively, until the values of T-^ and Tp, are within the 
desired accuracy for two consecutive iterations. 

The convective heat transfer coefficient for the flow of fluid in the channel is obtained 
by using the correlation for forced flow in a circular pipe of equivalent diameter [12]. 

h = 0.0155 Pr°-" 4^ (66) 

° ^ Deo 

The values of the Prandtl number and Thermal conductivity of the fluid are cal- 
culated at mean fluid temperature. Total heat transferred to the fluid is obtained by 
summing up heat gains by all the channels. 

For the case of high temperature working fluids, heat transfer coefficient are given 
in the form of. 


h = 1.62 


V, 


/ 


I D. 


f 


eq ^ 


Where / is given for different temperatures in [.30], 
flow inside tubes. 


( 67 : 


for the case of laminar 


4.6 Convective Heat Loss From The Plates 

For the ease of the analysis all the absorber plates and glass sheets are assumed at a 
uniform temperature at surface. Empirical correlations have been used for calculating 
the convective heat transfer coefficients. 
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4.6.1 Correlation for heat transfer coefficient from Glass sur- 
face to atmosphere 

For Convective heat loss from the outer side of the glass cover to the environment, 
ASHRAE Fundamental Volume [23] recommends the correlation for average coefficient 
of heat transfer for all orientation of surface 

= 5.50 -r2.70i;.> (68) 

4.6.2 Correlation for natural convection heat transfer coef- 
ficient from an inclined plate in rectangular enclosure 

For heat loss from the heated plate inclined at an angle 3 from the horizontal, due to 
natural convection, the most recent and reliable correlation for comparative collector 
design studies, is given by Hollands [16], 


Nu = 


+ 


1.0 + 1.44 



/ Gr Pr cos{i3) 
V 5830 


1708 \ 

Gr Pr cos ( 3) ) 

i 1 (+) 




sm^-^(1.8/3)1708 
Gr Pr CO si 3) 


(69) 


Where only positive values of the quantities in bracket with the {+) exponent are used 
and are taken as zero for negative values. The correlation is valid for the range of 8 
from 0 to 75 ° of inclinations. 
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4.6.3 Correlation for natural convection heat transfer coef- 
ficient from vertical heated and cooled plates 

The correlation given by Churchill and Ghu [27] is. 

Nul 

for 0 < Rai < 10^ 

for 10^ < Rcii 

.\11 properties except 3 in Or are to be calculated at mean film temperature i.e. 
^ T.-urfc.ce + Tftu.d ^ calculated at Tfi^id- 

4.6.4 Correlation for natural convection heat transfer coef- 
ficient from inclined heated and cooled plates 

Rich [21] recommended the use of g co.s{0) for plates inclined at an angle 9 from 
vertical. His results are substantially correct for lower surface of a heated plate or the 
upper surface of the cooled plate. Thus convective heat transfer coefficient for wall 2 is 
obtained using this correlation. For the upper surface of a heated plate and the lower 
surface of the cooled plate, the boundary layer becomes unstable and separates at low 
Gr values. This has been studied experimentally by Fuji! [25]. For this unstable case, 
he recommended the following correlations, 


= 0.680 + 0.670 i?aMl + 


0.492 \ ir 
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1 / 0.492 

0.825 + 0.387 Ral 1 + 
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1 


Nul 


0.13 Ra'l 

for 

Rai 

> 

5 10‘ 

.\ Ui 

= 

0.16 Ral 

for 

Rai 

< 

5 10 


<1) 


Where Gr is based on gcos{d). 


4.6.5 Correlation for natural convection heat transfer coef- 
ficient from vertical plate 

Chapman [21] recommends the following correlation. 


_ 0.637 Ra° -'^ 

A ^ Q 

1 ^ ~pr 

Where Lc is the ratio of .Area to Perimeter. 


(72) 


4.7 Modeling of Glass Sheet 

When a solar or thermal flux passes through the glass sheet, it is absorbed in the sheet 
according to Beer's law. Thus absorption is not uniform throughout the thickness. 
Similarly, the emission from the glass sheet is also a bulk phenomena, which causes 
self- absorption of the emitted radiation. Thus the analysis of the glass sheet is done 
by dividing the thickness in strips of equal thickness. Energy balance equation is 
written for all strips and then solved with the equations of the energy balance for other 
parts of the system. 
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4-7.1 Radiation Emitted By Glass Sheet and Strip-wise Ab- 
sorption of Primary Radiation 

The Beer's law has dependence on the direction of the radiation, because effective path 
length that the radiation travels is 


L p ““ 


Ax 


r3) 


cos{9) 

To account this direction dependence, absorption and emission from each strip is 
calculated in different directions and then summed up using .Simpson's rule to get the 
total absorption and emission. 

If the end effects are neglected, emission from a strip of thickness dx per unit solid 
angle du: is calculated using Equation 27 . 


3E = 


(T4) 


In the direction (9. the value of duj is. 


= sin(9)d6d6 (75) 

Where 9 varies from 0 to f and o from 0 to 2". 

.Since emission is not dependent on angle Equation 75 can be rewritten as, 

= 'lit sin{9)d9 ( 76 ) 

Using this definition in Equation 74 . 


(") 


dE = ‘2''itknlf^WBsin{d) 89 dxAi, 



CHAPTER 4. MATHEMATICAL MODELLING 


54 


If the thickness dx is of the order of the inverse of the absorption coefficient, it 
is fairly correct to assume that no self absorption occurs inside the strip. Thus for a 
strip of finite thickness, the emission from each face can be calculated by doing inte- 
gration over the strip thickness. Since the integral in 9 direction can not be calculated 
analytically. Simpson’s rule is here used to integrate in the 9 direction. 

Thus emission from each side of the strip of thickness Ax in direction 9 is. 


AEi{9) = 2~ithO’n\hAf,A9 [ THx)sin{9)€ ( 

Jo 


(7S\ 


Solving this, with the assumption that temperature inside the strip Tg{x} remains 
same throughout the entire volume of the strip, we get 


AEi{9) = 2an^t^AkA9THx)tan{9) 


-(zhAA] 

l—e 


179) 


Now, numerical integration in the direction 9 is done to calculate the total emission 
from the strip. Using Simpson’s rule, 


= - 


Tn—N0/2 


AEi{9 = 0)+AEi{9 = -)+2 ^ {AEi{{2m - l)A9) + 2AEA2mA9]} 

(80 


m=l 


Now this emitted radiation from strip is absorbed by the strips coming in its 
way to the surface. Thus the energ}’ absorbed in the strip from the energ}- emitted 
by all other strips is. 



52 AEi{9)e I ; 1-e 


( 81 ) 
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Where Dij is the center to centre distance between and strips less the thickness 
of one strip. 

For getting total amount of the absorption i.e. Simpson's rule is applied in 

the same manner as applied to AEi. This gives the absorption of the primary radiation 
emitted by \-arious strips only. For getting the emission that goes from glass to outer 
medium, reflectivity calculations at the glass air or glass-gas interface are required. The 
emission those which are incident at angle more than the critical angle are completely 
reflected back and through multiple reflections, get absorbed in the bulk itself. 

So firstly, the primary emission and from each strip reaching the interface 
is calculated and than the effect of interface reflectivity is accounted to get the total 
radiation emitted by the upper and lower surface of the glass sheet ( Qr Qi ) 
the absorption of the reflected primary radiation in each strip. 


Qf(e) = 

QriO) = 


X;AE.(9)e' 

f:^Ed6)e 


y cos{e) j 
I cos{ 6 ) J 


(82) 


(83) 


Where Di^u Is the distance between upper face of the strip and the upper surface of 
the sheet. Similarly. £),■,; is the distance between lower face of P^ strip and the lower 
surface of the sheet. 


4-7.2 Strip- wise Self-absorption of Secondary radiation 

Once the amount of the primary radiation reaching the interface in the direction 6. is 
known, multiple reflections of this is to be calculated. Transmissivity, reflectivity in 
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the direction 6 is calculated using Equation 24 . with a slight modification, that for 
angles greater than the critical angle reflectivity is unity. 


If 0 > then. 


sin 
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nth 


(84) 


li 9 < 6 ct then. 
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( 86 ) 
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frac-2^ - 

(1-p;^ 


/rac2„ = 


(i - 

frnci = 

O.o(/raci^ -f fraci ,, ) 

frac2 = 

0.5( /rac>_^ -f frac2 ., ) 


Q“(0) = frac,Q^;(0]~ frac.Qfiff) (88) 

QtiO) = fTaciQ^^(6) A f7'ac2Q‘^{d) 


In order to get total emission from each surface of the glass Simpson's rule is used to 
integrate in the 9 direction from 9 = 0 to 9 = 


q: 


m ,=.\' e /2 


Q' = o 


Q“(0 = O) + Q:1(^ = ;^) + 2 {Q“((2m-l)A0)+2(5“(2mA^)} 

“ m=:l 

__ m = Sgl 2 
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m=l 


For getting the energy absorbed in different strips due to secondary radiation i.e. 
multiply reflected part of the primary radiation from the interfaces, expression for 
energy absorbed in strip is derived on the similar lines. 

If 0 > then. 


31 


= 0.0 


^3. = 


( 39 ) 


0.0 
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1 - [piAG' 

Pan +-^311 '9 GuJ 


,1 - ('>;•, f-;") 

O.oi fraci^^ + /rac>,^ J 

0.5(/rac2^^ + /?'aci„ J 


(91) 


So. secondary energy absorbed in the strip is. 


Q"bs,(^) = Qr (^) + /^aC2J (6) 


( 92 ) 


The above equation is integrated ( using Simpson’s rule ) in the 0 direction to get 
the total secondary absorption in strip. 
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4.7.3 Strip-wise Absorption of Radiation Falling On Glass 
Sheet 

When the radiation passes through the glass sheet, it is absorbed in different amounts 
at different locations depending on the length of path traveled by the radiation to reach 
the location and absorption coefficient 7. 
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= 
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So. amount of beam solar radiation ( in direction 0 ) absorbed in strip is. 



1 - e 


L cos{ 6r ,) 



(96) 


This gives the absorption of the radiation in the direction 0. In the case of diffuse 
solar radiation the above quantity is numerically integrated over 0 = 0 to j. .Absorption 
of the thermal radiation is calculated using the same equations with jth instead of ''/so 
and Uf/i instead of rzso- 

It gives the distribution of absorbed radiation falling on the outer side of the glass 
sheet, with strip being the innermost and the outermost. 


distribution of beam solar radiation 
distribution of diffuse solar radiation 
0(1^ distribution of diffuse thermal radiation 
distribution of diffuse infrared radiation 


4.7.4 Energy Balance Equations For Glass Covers 

Once the absorption of the incident diffuse and beam solar radiation and thermal 
radiation in different strips, and emission of the thermal radiation from different strips 
and strip-wise self- absorption of the emitted thermal radiation are knowm. equations 
can be written to get the temperature increment of different strips after time At. 

For outer glass cover 

If there is only one glass covers, different radiations falling on the outer cover have 
been shown in Figure 11 . 

For lowermost strip i.e. j = 1, 
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Figure 11: Incident Radiations on Outer Glass Sheet When Single Cover 
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For uppermost strip i.e. j = Nx. 
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5.g=l,j = .Vi 






* 1 ,j=,Vr 

■\sb T^sb 


(ibsig-i^j^y^ 

th /^th 


A£',3=i._,=,Vx + QskyQ'd^^y^ (98) 
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For intermediate strips i.e. j = 2 to j = — 1, 
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Figure 12: Incident Radiation on Outer Glass Sheet When Multiple Covers 


NT 

Nt 


If there is more than one glass covers, different radiations falling on the outer cover 
have been shown in Figure 12 . 

For lowermost strip i.e. j = 1. 


SO 


= + QskyQ^: + 


(99) 
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For uppermost strip i.e. j = N^, 
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Figure 13: Incident Radiations on Inner Glass Sheet When Multiple Covers 
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For intermediate strips i.e. j = 2 to j = — 1. 


^’^3ig=l, j — [Qlbs,g=i,j + '54s.j=i,j — AAt-p=l,i + QskyQin, d- Q^n^ f^Qd; 
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[ 102 ) 
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ith 

K,A, 

Dx 


_J]/(M,Cp,(—)) 


For inn er glass cover 

Different radiations falling on the inner cover have been shown in Figure 13 . 


For lowermost strip i.e. j-= l~ 




CHAPTER 4. MATHEMATICAL MODELLIXG 


64 


AT 


3 i -V cover . j = i 
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- AE..= 
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(103) 


At 
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For uppermost strip i.e. j = .V^ 
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For intermediate strips i.e. j = 2 to j = A^j; — 1. 
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For intermediate glass covers 

Different radiations failing on the intermediate covers have been shown in Figure 


PL to 
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Figure 14: Incident Radiations on Intermediate Glass sheet When Multiple Covers 
For lowermost strip i.e. j = 1. 
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For uppermost strip i.e. j = AC, 
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4.7.5 Stability Condition 

In the above equations, choice of Ax and A/ are not independent of each other. For 
the stability of the solution, coefficient of the Tg ^^ , must be positive. 

For strip and Ni^" strip, the stability criteria requires : 


(Ax)^ ^ Kg AMs 
M ^ M,Cp, 

For intermediate strips, the stability criteria requires : 


lllO) 


( Ax)-^ KgAMg 

A( ^ - M,Cp, 


1111 ) 


4.8 Energy Balance Equations for Walls and Gas- 
mixture 

4.8.1 FOR BOX 

The temperature of the wails and the gas mixture in the box, after time interval At 
can be obtained by doing Energy balance of the heat fluxes absorbed and emitted. For 
the walls of the box, these equations are as follows 
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At 
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4.8.2 FOR gap 

Similarly, equations for increment in the temperatures of the walls and gas mixture in 
the gaps between the two glass covers can be calculated as 
For gap, 


At 


AT 


At 


+ aZiGt^,iG,,iAi,,TM - 

for i = 2. 3,4, 5 


ig^rntx 


igA-l,b 


= («£,„ + QS,„ + (p7 + o”:.*, )0S„ 

+ ^igj + ^ig,T + ^ig,L ^igX 

5 

f I ^so \ / /^iso • /"rso \ ^ so /"iSO 

- (s + + ^^9.0 ~ l.^i9.J^ig,J 

j=2 

+ Ql„ + Ql,„ - i< + NK-i + ''4‘.o> 

+ ^2 {^i9j^i3,j'^ig,j ~ Q cv^g,2)\l ^ig,mixC Pig^n 


j=2 


(113) 



CHAPTER 4. MATHEMATICAL MODELLIHG 


68 


4.9 Front Heat Loss 

The front heat loss consists of the following losses, 

[Ij- Convective heat loss from the glass surfaces to atmosphere, 

[2]- Radiative heat loss from the glass sheet to atmosphere, 

[Sj- Reflection of incident solar radiation. 

[4] - Transmission of the reflected solar radiation from the inside 

surface of the collector. 

[5] - Transmission of the thermal radiation from the inside surfaces of the 

collector. 

4.10 Efficiency 

The efficiency has been defined as the ratio of the total heat collected by the heat 
extracting or working fluid to the total incident beam and diffuse solar radiation on 
the outer glass sheet. 


V = 


Total Heat Transferred To Working Fluid 






Qt 


sd 

in 


(114) 



Chapter 5 


Solution Procedure 

Step [ 1] Read input data. 

Step [ 2j Initialize Temperature data. 

Step [ 3] Initialize emissivity, absorptivity of gas as zero. 

Step [ 4] Calculate Shape factors and Geometric mean beam lengths. 

Step [ 5] Calculate Transmissivity, reflectivity and absorptivity of the 
glass sheet for thermal, infrared, beam solar and diffuse solar 
radiation. 

Step [ 6] Calculate stripwise absorption of thermal. Infrared, beam solar 
and diffusesolar radiation in glass sheet. 

Step [ 7] Calculate properties of air at ambient temperature. 

Step [ S] Formulate matrix form of equation for the analysis of the glass 
sheets. 

Step [ 9] IF ( no gas in box and single cover) THEN 

Step [10] Calculate solar energy falling on different walls. 
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Step [11] Calculate solar irradiation on different walls of the box. 

GOTO Step [19]. 

Step [12] IF (no gas in box and gap) THEN 

Step [13] Calculate the solar irradiation on wails of different gaps. 

Step [14] Calculate solar irradiation on different walls of the box. 

GOTO Step [19]. 

Step [1-5] IF (gas in the box) THEN Calculate emissivity. absorptivity 
(for thermal and solar radiation both) of gas in box. 

Step [16] IF (gas in the gap) THEN Calculate emissivity .absorptivity 
(for thermal and solar radiation both) of gcis in gap. 

Step [17] Calculate solar energv- falling on different walls. 

Step [18] Calculate solar irradiation on different walls of the box and 
gaps. 

Step [19] Calculate heat transferred to the fluid. 

Step [20] Calculate heat radiated by the glass sheets. 

Step [21] Calculate secondary absorbed energy in glass sheets. 

Step [22] Calculate convective heat transfer from different walls of box 
and covers. 

Step [23] Calculate thermal irradiation on different walls of box and 
covers. 

Step [24] Calculate temperature rise in different strips of glriss sheets. 

Step [25] Calculate temperature rise in different walls of box and gaps. 
Step [26] Check time passed. IF (time < duration) GOTO Step [15] ELSE 
GOTO Step [27], 
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Experimental set-up 


6.1 Description 

Main components of the solar collector, fabricated in the lab for conducting experiment, 
are ; 

(1) Absorber plate 

(2) Flow passage 

(3) Glass Cover 

(4) Insulation 

(5) Containing box or casing 

(6) Gas mixture used in the box 

6.1.1 Absorber plate 

The absorber plate is fabricated from 24 gauge copper sheet. The choice of copper 
is based on high thermal conductivity to validate the assumption of uniform absorber 
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plate temperature. The absorber plates form a rectangular cavity with total five num- 
ber of absorber plates. The sheets are blackned on the face receiving solar radiation. 
Black oxide is sprayed on the sheets to increeise the solar absorptivity of the absorber 
plates. 

6.1.2 Flow passage 

The flow passage consists of rectangular channels of G.I. sheet, attatched to the side 
absorber plates in a spiral fashion, with two different flow paths. Selection of two 
flow paths is based on the requirement to e.xtract the maximum amount of heat and 
to keep the side absorber plate temperatures approximately same for all four plates. 
The channels are attatched to the base in the same manner, so that flow leaves the 
absorber plate at the centre of the base absorber plate. This is done to extract the 
maximum amount of heat, because the centre is expected to be free from corner effects 
and to have the maximum temperature. Thus, it has been tried to keep the difference 
between the absorber plate and fluid temperature minimum and the fluid should come 
into contact with the wall at progressively higher temperatures. 

6.1.3 Glass Cover 

The experimental set-up is fitted with two number of glass sheets of thickness 4mm. 
The thickness is chosen as a compromise between the strength of the sheet and the 
heat loss. The gap between the two glass sheets is kept 1 inch. As the gap thickness 
is increased, the total heat losses are increased, hence the gap is kept at the minimum 
but with keeping the fabrication difficulties in mind. Rubber beading and packing 
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rubber sheet are used to make the gap air tight. Properties of different possible cover 
materials are given in all Solar Energy Handbooks [13]. 

6.1.4 Insulation 

Glass wool insulation is used to insulate the flow channels in order to reduce the heat 
losses from side and base through the casing walls. Glass wool is selected because it 
can easily withstand temperatures of the order of 200 — 300 °C. 


6.1.5 Containing box or casing 

The absorber plates fitted with flow channels and insulated with Glass wool is kept in 
the G.I. casing of rectangular shape. The selection of the G.I. sheet is based on the 
strength and heat loss considerations. 

6.1.6 Gas mixture used in the box 

Ammonia gas is filled in the box. For filling and evacuating purpose two passage have 
been given. The lower passage is used for filling purpose whereas upper passage is 
used for evacuating. During filling, liquid ammonia is put in the box through the lower 
passage and upper passage is kept lock, air being heavier than ammonia, leaves out 
from the lower passage. The exact composition of the gas is difficult to measured, but 
an approximate can be made easily. To evacuate the gas mixture from the box. both 
passages are kept open, resulting in the leaving out of the gas-mixture through the 
upper passage and entering of fresh air from the lower passage. 
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6.1.7 Heat extracting fluid used 

Depending upon the range of temperatures in operation, different fluids can be used 
for this purpose. Heat transfer fluids along with their density, specific heat and boiling 
point, heat transfer coefficient parameters are given in it30:. For the temperatures below 
100 °C. water is the best suited collector fluid. 

6.2 Solar radiation Intensitv Measurement 

%/ 

Solar radiation intensity on horizontal surface was measured using a Kipps and Zonen 
CM5 pyranometer (See Appendix C. Figure 36 ). It was put near to the collector 

and leads were connected to the digital milivoltmeter. Thermo-emf generated was 
measured and converted into solar inntensity by multiplying it by the pyranometer 


constant. 
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Results and Conclusions 


7.1 INPUT DATA for Computation 

[1]. Geometric parameters : 

1,= Liy= 60.0 cm. 

I.-= 15.0 cm. 

A 5 = 4.0 mm. 

■^gap= 3.0 cm. 

Lcha= 3-0 cm. 

3= 40.0 ° 

[2]. Optical properties of glass ; 

"th— 8000.0 

4400.0 

lso= 30.0 m~^ 
nth- 1-35 
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nin= 1.40 
nso= 1-52 

[3] . Thermal properties of wall material : 

Kcopper= 386.0 W/rn — A' 

386.02 J/Kg. - K 
3/copper = 4.613 Kg-lm.- 

^painiedcopper — 0.9*5 

(y^° ■ = 0 98 

paint edcopper ^ 

_ A qrr 

^ paintedcopper 

[4] . Thermal properties of Glass : 

I'^'giass= 0.80 W/rn - K 
Cp,u....= 669.92 J/Kg. - K 
Tigiass= 10.000 Kg./m^ 

[5] . Average Thermal properties of air : 

<^Pa.r= 1004.88 J/Kg. - A' 

Al — 1.1650 Kg./m^ 

[6] . Average Thermal properties of A'As : 

a = 2m.60 J/Kg. - K 
M = 0.6880 Kg./m^ 

[7] . Solar parameters : 

A<ioy= 115 
Lattitude= 26.6 ° 
hc= 1353.0 W/m- 
Kt= 0.80 
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a(Kx)= 0.270 
6(AV)= 0.395 

[8] . Ambient air parameters : 

VaiT= 3.00 m/s 
RH= 0.50 
CC= 0.00 
a(sky)= 0.6164 

[9] . Modelling parameters ; 

Ab= 16 
.V5= 180 
At= 0.05 

TimePeriod— 2 Hours ( 1ST 11:00 to 13:00 Hours ) 

7.2 Computational Results for Nontracking col- 
lector 

The computational results for the nontracking collector with the parameters listed in 
the above section, can be explained in two sections, first being the transient perfor- 
mance of the collector and the second the overall performance of the collector. 

7.2.1 Transient Performance of Nontracking Collectors 

For the study of the effect of the absorbing media on the transient performance of the 
collector, results have been computed for partial pressure of NH^ = 0.0, 0.4, 0.8. 1.0 
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time in minutes 

Figure 15: Xontracking ; Variation of Input Solar Energy with time 

( on the base plate ) 

and evacuated case with no selective coating. The results have been computed for the 
fixed geometry of the collector and water flow rate. 

Ln = Lh = 60.0 cm. 

Lc = 15.0 cm. 

Lgap = 3.0 cm. 

Vw = 0.005 m/s 

7.2. 1.1 Variation of Input Solar Energy with time 

Input solar energy consists of two parts, diffuse and beam solar radiation. The beam 
component is proportional to the cosine of the incidence angle of the falling radiation. 
The incidence angle varies with the hour angle according to Equation 3 . The diffuse 
component depends upon the hour angle because of the hour angle dependence of the 
extra-terrestrial horizontal solar intensity. The variation of input solar energy with 
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Figure 16: Nontracking : Variation of Working Fluid Outlet Temperature with time 
time is shown in Figure 15 . 

7. 2. 1.2 Variation of Working Fluid Outlet Temperature with time 

The angle of incidence of the beam radiation changes from zero at solar noon to finite 
values for before and after solar noon. Thus the radiation falling on the base plate 
increases before the solar noon and starts decreasing after the solar noon, which results 
in the fluid outlet temperature rise before solar noon and fall after solar noon. The 
variation of fluid outlet temperature with time is shown in Figure 16 . 

7. 2. 1.3 Variation of Outer Glass Temperature with time 

The variation of the outer glass temperature depends on the solar and thermal radiation 
falling on it and optical properties of the glass. Since the use of the NH^ reduces the 
thermal flux on the glass sheet, its temperature is reduced with increasing partial 
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Figure 17: Nontracking : Variation of Outer Glass Temperature with time 
pressure. The variation of Outer glass temperature with time is shown in Figure 17 . 

7. 2. 1.4 Variation of Front Heat Loss with time 

Front heat loss consists of two components, first is the radiative loss from the outer 
side of the glass sheet and second is the convective heat loss from the outer side of 
the glass sheet. The convective heat loss coefficient strongly depends upon the wind 
velocity. The variation of front heat loss with time is shown in Figure 18 . 

7.2. 1.5 Variation of Heat Extracted by Working Fluid with time 

The amount of heat collected is directly proportional to the incoming energy. Thus 
the trend of this curve is similar to that of the total input solar energy. The variation 
of heat extracted by working fluid with time is shown in Figure 19 . 
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Figure IS: Nontracking : Variation of Front Heat Loss with time 
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Figure 19: Nontracking : Variation of Heat Extracted by V^orking Fluid with time 
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Figure 20: Nontracking: Collection Efficiency vs Fluid outlet temperature 

7.2.2 Overall Performance of Nontracking Collectors 

For the study of the effect of the absorbing media on the overall performance of the 
collector, results have been computed for partial pressure of .V/f 3 = 0.0. 0.4. 0.8. 1.0 
and evacuated case with no selective coating. The results have been computed for the 
fixed geometry of the collector and different water flow rates. 

- 7.2.2.1 Variation of Collection Efficiency with Fluid outlet temperature 

The plot of the p vs Figure 20 , for various partial pressure of the gas in the 

collector box, clearly indicates that the use of the absorbing media i.e. N H 3 increases 
the collection efficiency by fairly good percentage as compared to the air filled collector 
at high temperatures. At lower temperature the temperature of the absorber plate 
(base) is not very high, giving the case of enclosure with uniform wall temperatures. 
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Figure 21: Nontracking : Total Front Heat Loss vs Partial pressure of N Hp, 

Thus radiation flux falling on all sides of the box are nearly same, resulting in marginal 
gain in the efficiency. The use of the absorbing media cuts down the high intensity 
radiation flux on the glass sheet from the inner side ( due to higher temperature of 
the base plate ) to the glass sheet by absorbing it and then re-emits it uniformly in all 
directions, causing a reduction in the thermal flux on the glass sheet from inner side. 

The variation of the increase in the efficiency with partial pressure depends upon 
the emissivity and the absorptivity of the gas. It increases with partial pressure of 
N Hz. at faster rate in the starting and becomes almost constant at higher values. 

7. 2. 2. 2 Variation of Total Front Heat Loss with partial pressure of NHz 

The total front heat loss decreases with increasing partial pressure. But the difference 
in the values becomes smaller and smaller with increasing flow rates, i.e. at lower fluid 
outlet temperatures. This happens because at lower absorber plate temperature both 
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Figure 22: Nontracking : Total Heat Extracted by Working Fluid vs Partial pressure 
of /Vf/3 

the glass sheet and the absorber plate is at nearly same temperature, which makes the 
use of gas ineffective. The variation of total front heat loss with partial pressure of 
N is shown in Figure 21 . 


7.2. 2.3 Variation of Total Heat Extracted by Working Fluid with partial 
■ pressure of NH^ 

The total heat extracted by working fluid increases with increasing partial pressure. 
But the difference in the values becomes smaller and smaller with increasing flow rates, 
because of the reason stated in the case of the total front heat loss. The variation of 
total heat extracted by working fluid with partial pressure of NH 3 is shown in Figure 
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Figure 23: Nontrar.king : Outer Cover Temperature vs Partial pressure of N Hn 


7. 2. 2.4 Variation of Outer Cover Temperature with partial pressure of .V //3 


The outer cover temperature increases with increasing partial pressure. But the differ- 
ence in the values becomes smaller and smaller with increasing flow rates, because of 
the reason stated in the case of the total front heat loss. The variation of outer cover 
temperature with partial pressure of N Hz is shown in Figure 23 . 


7.2. 2 . 5 Variation of Working Fluid Outlet Temperature with partial pres- 
sure of NHz 

The variation of working fluid outlet temperature with partial pressure of NHz is shown 
in Figure 24 . It also shows the same behavior that, at lower temperatures there is 
only marginal gain in the outlet temperature whereas at higher outlet temperatures 
gain is sufficient. 




CHAPTER 7. RESULTS AND CONCLUSIONS 


86 



Figure 24: Nontracking : Working Fluid Outlet Temperature vs Partial pressure of 
NHs 

7.3 Computational Results for Tracking Collector 

The computational results for the tracking collector with the parameters listed in the 
above section are explained as follows. The incident beam radiation is constant in this 
case because of the tracking of the collector. Hence all quantities attain a steady state 
value after some period. Thus in the transient study values first increase rapidly and 
with the passage of the time they achieve steady state value. 


7.3.1 Transient Performance of Tracking Collector 

For the study of the effect of the absorbing media on the transient performance of the 
collector, results haveTieen computed for partial pressure of NHz = 0.0. 0.4, 0.8. 1.0 
and evacuated case with no selective coating. The results have been computed for the 
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Figure 25: Tracking : Variation of Time for achieving steady state with Partial pressure 
of V //3 

fixed geometry of the collector and water flow rate. 

7. 3. 1.1 Variation of Time for achieving steady state with Partial pressure 
of NH 3 

The variation of time required to attain the steady state by the collector with partial 
pressure of N Hz is shown in Figure 25 . Time taken by the collector to achieve the 
steady state is least in the case of the evacuation and it increases with the partial 
pressure of the N Hz, because the gas itself starts taking acti^^e part in energy transfer. 

The values of the front heat loss, heat transferred to the working fluid, working 
fluid outlet temperature and the outer glass temperature at the steady state had been 
plotted and they show the same trend as in the case of nontracking collectors. 
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Figure Tracking ; Collection Efficiency vs Fluid Outlet Temperature 

7.3.2 Overall Performance Of Tracking Collector at steady 
state 

For the study of the effect of the absorbing media on the overall performance of the 
collector, results have been computed for partial pressure of N Hs = 0.0. 0.4. 0.8. 1.0 
and evacuated case with no selective coating. The results have been computed for the 
fixed geometry of the collector and different water flow rates. 

The plot of the rj vs Tji^,. Figure 26 , for various partial pressure of the gas in the 
collector box, clearly indicates that the use of the absorbing media i.e. N Hz increases 
the collection efficiency by fairly good percentage as compared to the air filled collector 
at high temperatures. At lower temperature the temperature of the absorber plate 
(base) is not very high, giving the case of enclosure with uniform wall temperatures. 
Thus radiation flux falling on all sides of the box are nearly same, resulting in marginal 
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figure 27: Tracking : Total Heat Extracted by Working Fluid vs Partial pressure of 
V//3 

;ain in the efficiency. The use of the absorbing media cuts down the high intensity 
adiation flux on the gleiss sheet from the inner side ( due to higher temperature of 
he base plate ) to the glass sheet by absorbing it and then re-emits it uniformly in all 
irections. causing a reduction in the thermal flux on the glass sheet from inner side. 

The variation of total heat extracted by working fluid with partial pressure of .V //3 
> shown in Figure 27 . 

The variation of total front heat loss with partial pressure of N H:i is shown in 
‘igure 28 . 

The variation of outer cover temperature with partial pressure of N Hz is shown in 
igure 29 . 

The variation of working fluid outlet temperature with partial pressure of NHz is 
lown in Figure 30 . 
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igure 28: Tracking : Total Front Heat Loss vs Partial pressure of N H 
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Figure 29: Tracking : Outer Cover Temperature vs Partial pressure of NH 
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Figure 30: Tracking :Working Fluid Outlet Temperature vs Partial pressure of .V //3 

The explanation of the above plots is the same as that given in the case of non- 
tracking collector. 


7.4 Experimental Results 

Experiment? had been conducted for the polar tracking mode with gap between the 
cover and the absorber plate filled with air and gas respectively. The values of the 
solar intensity, ambient temperature and fluid inlet and outlet temperature and flow 
rate of the working fluid had been measured. Finally the collection efficiency had 
been calculated and plotted. Figure 31 , with outlet temperature. The plot show 
the same behavior as predicted by the computational results, but the values are not 
matching exactly. The main reasons for the difference in the values lies with the 
optical properties of the glass, convective heat transfer in the channel and variation in 
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igure 31: Experimental Results: Collection Efficiency vs Fluid Outlet Temperature 


3 ambient conditions. 


5 Effect of Gap between Glass Sheets on Collec- 
tor Performance 

the distance between the two glass sheets is increased, less thermal radiation reaches 
the outer glass sheet, because the side walls also starts coming into picture bv 
eiving some amount of the radiation emitted by the lower glass sheet. At the same 
le in the case of the nontracking collector, the side walls starts intercepting the solar 
liation, thereby, reducing the total solar energy input to the absorber plate. This 
i be avoided if the channels are attatched to the side walls starting from the top, i.e. 
ter also extract heat from the space between the two glass sheets. In the present 
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Figure 32: Nontracking : Working Fluid Outlet Temperature vs Gap 


case the analysis as well as experimental set up uses this configuration. If the gas is 
used in the gap also, the increase in the gap causes more absorption of the thermal 
radiation and reduced emission towards the outer glass sheet because of the reduced 
configuration factor. 


These factors result in the overall increase in the efficiency and the outlet tempera- 
ture of the working fluid. The increase in the efficiency is increasing with gap distance 
initially and starts decreasing when the gap distance is too much. The computational 
results have been plotted for both the cases, nontracking, Figure 32 , 34 , and tracking, 


Figure 33 , 35 . 
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Figure 33: Tracking : Working Fluid outlet Temperature vs Gap 



Figure 34: Nontracking : Heat Collection Efficiency vs Gap 
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Figure 35: Tracking : Heat Collection Efficiency vs Gap 

7.6 Conclusions 

From the computational and experimental results, the following conclusions can be 
made : 

[1] Use of the gases which are transparent for solar radiation and strongly absorb- 
ing and emitting for thermal radiation, can be helpful in improving the collector's 
performance. 

[2] But their use is effective only at higher temperature applications. .A.t lower 
temperatures, the improvement is not sufficient to be of practical importance. 

[3] Use of the absorbing media can be helpful in improving the performance of 
the concentrating type of solar collectors, where absorber plate temperature is usually 
much higher than the glass cover temperature. 

[4] Use of the absorbing media can also be helpful in improving the performance of 
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the nonconcentrating type of solar collectors, where reflector mirror is used to increase 
the intensity of solar radiation to achieve higher temperatures. 


Table : Readings for collector with air in it : 


'low rate 

Tambient 

Heat input 

Tfluidint^, 

Tfluid^ml^f 

Heat Extracted 

Efficiency 

1 1pm 

"C 

Watts 

°c 

°C 

W atts 

% 

.167 

42.0 

325.56 

32.0 

49.0 

197.71 

60.70 

.100 

42.0 

325.56 

32.0 

56.0 

167.47 

51.40 

.059 

42.0 

325.56 

32.0 

63.0 

135.46 

41.60 

1.031 

42.0 

325.56 

32.0 

78.0 

100.33 

30.80 

Table : Readings for collector with gas in it : 

^low rate 

Tambient 

Heat input 

TfluidinUt 

Tfluidoutlct 

Heat Extracted 

Efficiency 

n 1pm 

°c 

Watts 

°c 


Watts 

% 

3.222 

41.0 

319.30 

33.0 

46.0 

217.10 

68.00 

3.117 

41.0 

319.30 

33.0 

54.0 

172.40 

54.00 

3.074 

41.0 

319.30 

33.0 

62.0 

149.90 

47.00 

0.036 

41.0 

319.30 

33.0 

76.0 

112.15 

35.10 



Chapter 8 


Suggestion for Further Study 


• T 1 k‘ slucly can tx' oxteiided to concentrating type of collectors, where it is ex- 
p('ct-('d 1.0 giv(' ('ven better results. 

• Tlu' st udy can b(> ('xtended for other absorbing-emitting gases which are trans- 
parcmt for solar radiation and highly absorbing for thermal radiation. 

• The study can be extended to include the effect of spectral emission of the gaseous 
media, because in the case of selective surface use of the Hottel’s emissivity chart 
ar<' not valid. Ev('n for present case it is not very accurate. 

• l’h(' absorlx'r plat(‘ temperature had been assumed uniform throughout the plate, 
which is not. strictly true. The effect of this can be studied by dividing the plate 
in number of zones. It will increase the computational complexity, but will gi've 
more accurate results. 
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Tables Related to Solar Intensity Calculations [ 28 ]. 
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Appendix: A 


0 .') 


Table <): ( 'oeffi.'ients n{ hr) and /v» for Beam Radiation Model 

Interval 

aiky) 


0.00, 0.05 0.04 0.00 

0.05. 0.!.^ 0.01 0.002 

0.15,0.25 0.06 -<■‘006 

0.25. 0.35 0.32 ^0.071 

0.35. 0.45 0.82 -0.246 

0.45.0.55 1.56 -O.5-9 

0.55, 0.65 1.69 -0.651 

0.65.0.75 1.49 -0.521 

0.75. 0.85 0.2: 0.395 


I'ahU' 10: Coeffirient t/,;^ lor Sky Flux Equation 
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pendix: B 

p,elated to Geometric Mean Beam Lengths [10). 
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ible 11: Geometric Mean Beam Length Ratio's and Configuration Factor for Parailei 
Equal Rectangular Plates 


rjsaale, f** Configuration Faaor, r»«Gcornetric Mean 
Beam Length 




0 

O.I 

0.2 

0.4 

0.6 

0 

rie 

f 

1.000 

I.OOl 

1.003 

1.012 

1.025 

0.1 

r/c 

1.001 

1.002 

1.004 

1.013 

1.026 


F 


0.00316 

0.00626 

0.01207 

0.01715 

0.2 

r/c 

1.003 

1.004 

1.006 

1.015 

1.028 


F 


0.00626 

0.01240 

0.02391 

0.03398 

a4 

r/c 

1.012 

1.013 

1.015 

1.024 

1.037 


F 


0.01207 

0.02392 

0.04414 

0.06560 

0.6 

r/c 

1.025 

1.026 

1. 028 

1.037 

1.050 


F 


0.01715 

0.03398 

0.06560 

0.09336 

1.0 

r/c 

1.055 

1.056 

1.058 

1.067 

1.080 


F 


0.02492 

0.04941 

0.09554 

0.13627 

2.0 

r/c 

1.116 

1.117 

1.120 

1.129 

1.143 


F 


0.03514 

0.06971 

0.13513 

0.19342 

4.0 

r(e 

1.178 

1.179 

1.182 

1.192 

1.206 


F 


0.04210 

0.08353 

0.16219 

0.23271 

6.0 

rfc 

1.205 

1.207 

1.210 

1.220 

1.235 


F 


0.04463 

0.08859 

0.17209 

0.24712 

1 0.0 

r/c 

1.230 

1.232 

1.235 

1.245 

1.261 


F 


0.04671 

0.09270 

0.18021 

0.25896 

20.0 

r/c 

1.251 

1.253 

1.256 

1067 

12182 


F 


0.04829 

0.09586 

0.18638 

0.26795 

03 

r/c 

Ijr72 

1.274 

1.277 

1.289 

U06 


F 


0.04988 

0,09902 

0.19258 

0.27698 
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Kipps and Zonen CM5 Pyranometer 
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Appendix: C 
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SOLARIMETER 

FOR OUTDOOR INSTALLATION 
CM 5 . CM 6 


The solarimeter base is provided with a spirit level. It 
should be placed in a horizontal position by means of the 
two levelling screws. 

The solarimeter CM 5 is screwed onto the base by means 
of the three screws provided. 

The screen is fixed by pushing the three studs into the 
three supports on top of the base and tightening the srews. 

The instrument should be positioned in such a way that the 
output cable is pointing NORTH. 

In order to avoid condensation on the inside of the glass 
domes the interior of the solarimeter is kept dry by means 
of a built-in drying cartridge. 

This drying cartridge (1) can be withdrawn from the mount- 
ing after removal of the retaining screw-ring (2) as indicated 
in fig. 1. Upon pulling the two parts of the cartridge apart, 
the perforated -tube can be filled with new silicagel or other 
suitable drying agent. 

For cleaning, the outer glass dome can be taken off as 
indicated in fig. 1. Upon unscrewing the retaining screw- 
ring (3) the glass dome (4) can simply be lifted off. In re- 
assembling care should be taken that the ring is screwed 
home very tightly. 


Connection to measuring equipment 

The solarimeter is provided with a two-core output cable. 
Black is the negative terminal and blue the positive. 

If long cables are used between the solarimeter and the 
measuring equipment, the possible influence of the resist- 
ance of these cables should be taken into account. 

In the case of measuring equipment with a high input 
resistance (such as the integrator CC 1) the resistance of 
the connecting cable can usually be neglected. 

When recording meters with a low internal resistance are 
used, the risistance of the cables may have a considerable 
influence. A correction is therefore required in order to 
avoid faulty measuring results. 

For the calibration in conjunction with the solarimeter of an 
indicating instrument not supplied by us. two methods can 
be applied as shown in figures 2 and 3 page 4. 

These schemes also enable the operator to determine the 
eventual influence of the resistance of the connecting 
cables. 



Figure 36; Kipps and Zonen CMS Pyranometer 




